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We study the interplay of magnetic and superconducting order in single crystalline hole doped 
Bai_a;Naa;Fe2As2 using muon spin relaxation. We find microscopic coexistence of magnetic order 
and superconductivity. In a strongly underdoped specimen the two forms of order coexist with- 
out any measurable reduction of the ordered magnetic moment by superconductivity, while in a 
nearly optimally doped sample the ordered magnetic moment is strongly suppressed below the su- 
perconducting transition temperature. This coupling can be well described within the framework of 
an effective two-band model incorporating inter- and intra-band interactions. In optimally doped 
Bai_a;Naa;Fe2As2 we observe no traces of static or dynamic magnetism and the temperature depen- 
dence of the superfluid density is consistent with two s-wave gaps without nodes. 

PACS numbers: 74.70.Xa, 61.05.C-, 74.62.Dli, 76.75.+i 



Strongly correlated electron systems like the cuprate 
high- Tc, organic, and heavy-ferniion superconductors ex- 
hibit a delicate interplay of superconductivity and an- 
tiferromagnetic order. Heavy-fermion superconductivity 
often takes place in the vicinity of a magnetic quantum 
critical point [1]. Similarly, magnetic order lies nearby 
and competes with superconductivity in the phase dia- 
grams of the cuprate high- Tc superconductors [2-4]. In 
the fcrropnictide superconductors the interplay of mag- 
netism and superconductivity leads to a new phase that 
is both homogeneously superconducting and magneti- 
cally ordered [5-8] — i.e., magnetic order and supercon- 
ductivity coexist microscopically [6, 9, 10]. Still, both 
instabilities compete for the same electron states around 
the Fermi surface, which causes a reduction of the or- 
dered magnetic moment upon cooling below Tc- In 
some materials electronic phase separation appears in- 
stead of microscopic coexistence. Experiments with elec- 
tron doped BaFe2-a;Coa;As2 suggest microscopic coexis- 
tence [9, 11, 12], nevertheless a recent report by Bern- 
hard et al. [13] indicates electronic inhomogeneity arc 
important near optimal doping. Electronic phase separa- 
tion, on the other hand, commonly occurs in hole doped 
Bai_2:K2;Fe2As2 [11, 14-16]. Only recently Wiescn- 
mayer et al. [17] showed a reduction of the ordered mo- 
ment by superconductivity in polycrystalline underdoped 
Bai_a;Ka;Fe2As2 — an unambiguous proof of microscopic 
coexistence — consistent with other reports [18-20]. 

In this Letter we report a strong interplay of the or- 
der parameters of superconductivity and magnetic or- 
der in Bai_2:Naa;Fe2As2. In a nearly optimally doped 
compound, the ordered magnetic moment of the Fe spin 
density wave order, that appears below Tn ~ 40 K, is 



reduced by k, 65% upon cooling below Tc ~ 30 K. 

We examined platclike Bai_^Naj;Fc2As2 single crys- 
tals (growth conditions and characterization arc de- 
scribed in Ref. [21]) with the c-axis oriented along the 
muon beam. The muon-spin relaxation was measured 
in zero (ZF) and transverse (TF) magnetic fields of 
up to 0.64 T using the DOLLY and GPS instruments 
at the Paul Scherrer Institute equipped with ^He flow- 
cryostats. To analyze the /iSR data, we used the 
MUSRFIT framework [22]. The stoichiometry of each 
crystal was determined with electron-dispersive x-ray 
spectroscopy (EDX). Magnetization, electrical resistiv- 
ity, angle-resolved photoemission spectroscopy (ARPES), 
Hall-effect, and specific- heat measurements of single crys- 
tals from the same batches are published in Refs. [21, 23, 
24]. 

Fig. la shows the muon spin polarization P{t) in ZF 
perpendicular to the beam. For x = 0.40, the relaxation 
of P{t) is temperature independent and characteristic 
for a paramagnetic material (down to the lowest mea- 
sured temperature of 2 K). For x = 0.25, and 0.30 we 
find two precession signals a ^ A,B with frequencies fa 
that indicate two magnetically nonequivalent muon sites, 
and prove long range magnetic order with T/v = 123(1), 
and 40.0(5) K, respectively. We analyze P{t) with two 
damped-cosine oscillations, and a nonrelaxing signal due 
to paramagnetic and non-oscillating signals. 

The muon-spin precession frequencies (which are pro- 
portional to the magnetic order parameter) of x = 0.25, 
and 0.30 are shown in Figs, lb, and Ic. For x = 0.25 we 
find no anomaly faiT) at Tc ~ 9 K on cooling (Tc from 
Ref. [21]). A fit with a general order parameter temper- 
ature dependence /a(T) = /a(0)(l - [T/TnTY, shown 
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FIG. 1. ZF /^SR of Bai-a:Nai:Fe2As2: (a) the muon spin polarization time dependence for x = 0.25 (left), 0.30 (middle), and 
0.40 (right); the temperature dependence of the magnetic order parameter (oc fa{T)) for (b) x — 0.25 does not show any 
anomaly at Tc, the inset shows the low temperature regime in which superconductivity occurred (the open symbols indicate 
the data excluded from the fit; the lines are fits, see text), but for (c) x = 0.30, the precession frequencies (the magnetic order 
parameter) are reduced by « 65% below superconducting critical temperature Tc2 (the line is calculated from an itinerant 
two-band model, see text) — nevertheless, all evidence supports microscopic coexistence in both cases. 



in Fig. lb, to the data above Tc reproduces the data also 
below Tc- By setting a = 1 we estimate /3 = 0.13(1), con- 
sistent with the 2D-Ising universality class (/3 = 0.125), 
and similar to results for BaFe2As2 [25]. 

We conducted fj.SK "pinning experiments" to verify 
bulk superconductivity: We cooled the sample with 
X = 0.25 to T = 5.0 K in a transverse magnetic field 
fJoHi = 500 mT, recorded a time histogram (Fourier 
transform shown in Fig. 2a as blue line), then isother- 
mally reduced the field to ij,qH2 = 450 mT and recorded 
another histogram (red line). The magnetic order caused 
four precession signals in the external magnetic field: 
7^/ioi?i ± 27r/^, and ± Stt/b, as indicated by 

the Fourier transform in Fig. 2a, in addition to the back- 
ground signal due to muons stopped by the cryostat walls 
or the detector. The signal fraction that followed the 
change of the external magnetic field was « 10% of the 
total signal, which is consistent with a background sig- 
nal. The remainder of the signal maintained its average 
internal field after we reduced the field, which indicates 



strong flux pinning due to superconductivity in the ma- 
jority of the sample. Instead, caused by enhanced vor- 
tex lattice disorder, only the damping of the spin preces- 
sion increased (spectral lines broadened). This evidence 
for bulk superconductivity taken together with the bulk 
magnetic order supports microscopic coexistence of mag- 
netic order and superconductivity, albeit with a coupling 
of the order parameters that is too small to be detectable 
[26]. 

The pinning experiment for x = 0.30, shown in Fig. 2b, 
also indicates bulk superconductivity. In contrast to 
X = 0.25, the magnetic order parameter for x = 0.30 
(see Fig. Ic) was suppressed by w 65% on cooling below 
Tc2 — 29.0(6) K, which proves microscopic coexistence 
of magnetic order and superconductivity and strong cou- 
pling of both order parameters. 

Theoretical studies show that microscopic coexistence 
as well as phase separation arise in itinerant multi-band 
models as a function of Fermi surface topology, band dis- 
persion, and symmetry of the superconducting order pa- 
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FIG. 2. Bulk superconductivity is indicated by TF-/^SR 
pinning experiments for (a) x = 0.25: After isothermally 
changing the field only ~ 10% of the signal precesses with 
/m ~ 7^/(21") ■ 450 mT due to muons that stop in the cryostat 
walls and the sample holder (bgd.), ±/q indicate signals from 
muon sites a = A, B with the local field parallel (+)/anti- 
parallel ( — ) to the applied field, and for (b) x = 0.30 the 
background signal amounts to ~ 5% of the total signal — both 
experiments indicate strong fiux pinning due to superconduc- 
tivity in the majority of the sample volume. 



rameter [5-8]. Within these models, s± Cooper-pairing, 
in which the gap changes sign between the hole and elec- 
tron pockets, has the largest tendency towards micro- 
scopic coexistence. In order to show that such an ap- 
proach can account for our measurements, we calculate 
the temperature dependence of the ordered magnetic mo- 
ment M(T) within a simplified two-band model incor- 
porating inter- and intra-band interactions that drive 
the spin density wave order (SDW) as well as super- 
conductivity (SC) [27, 28]: H = Eka{ec(k)cLck<T + 

e/(k)/L/k<T} + Ek,k',q,^.<T'{"l4+qa/k'-q<T'/kV'Ck^ + 

^(/k-Hq<T/k'-q<T'Ck'fT'Ck<T + H.c.)}, with the dispersion 
ec(k) of the electron (c) and ey(k) of the hole (/) band. 
Decoupling the interactions in mean-field we diagonal- 
izc the Hamiltonian by sequential Bogolyubov transfor- 



mations that lead to the eigenenergies Vl^ = [{E'^) -t- 
(A^) where 7 = a, /3 and E^'^ = e+±[(e-)'+M2]0-5. 
We use the following parameters appropriate for a two- 
band model with hole doping: t = 1.0, t' ~ 0.72, ef = 
1.3, Ec = -3.3, fi = 0.1, usDW = ui +U3 = 4.985, Usc = 
U3 — 5.605. Here M is the antifcrromagnctic (SDW) or- 
der parameter and are the two SC gaps (see Rcf. [28] 
for details). We solve the self-consistency equations in 
order to obtain the temperature dependence of the mag- 
netic order parameter M(T), and the gaps Ag(T), A^(r) 
of s± superconductivity. 

The results for M{T) and the superconducting gaps are 
shown as a solid/dashed lines in Fig. Ic. The dip in M (T) 
coincides with the onset of SC at Tc- In particular, the 
relative reduction of the magnetic order parameter below 
Tc is mostly determined by the magnitude of the super- 
conducting order parameter that competes for the same 
Fermi surface points. More elliptical electron pockets re- 
duce the magnitude of M[T) in comparison to the super- 
conducting order parameter, thus increasing the relative 
suppression. The agreement between the calculated tem- 
perature dependence and the experimental data is very 
good for T > 25 K. At lower temperatures the experi- 
mental data drop somewhat faster than the theoretical 
temperature dependence. The discrepancy may be con- 
nected to a change of magnetic fluctuations that change 
the magnetic moment at the SC Tc due to gapping of 
particle-hole excitations. 

For X — 0.25, the temperature independent total 
oscillation amplitude and TF-^SR show that 100% of 
the sample is magnetically ordered (not shown). For 
X = 0.30, the remaining precession amplitude in 5 mT 
TF /xSR indicates a paramagnetic phase with a volume 
fraction of 23(2)% below Tm, shown in Fig. 3a. The dia- 
magnetic shift of the average internal field of the para- 
magnetic phase already for T > Tc2 = 29.0(6) K in- 
dicates that this phase becomes superconducting below 
Tci = 34.5(13) K. The precession amplitude in zero field 
(see Fig. 3b) is temperature independent, which indicates 
that the magnetic volume fraction reaches « 77(2)% 
and remains temperature independent below w 39 K. 
The sharp increase of the magnetic volume fraction in 
Fig. 3b shows that, even though the transition tempera- 
ture is reduced to T/v = 40.0(5) K, the transition remains 
sharp — much sharper than in other ferropnictides (mea- 
sured from onset to saturation of the magnetic volume 
fraction it is not broader than 5(1) K) [13, 17, 29, 30]. 
The small damping rate At of the muon spin preces- 
sion also indicates well defined magnetic order: it did 
not exceed Xt/ {^irfA) ~ 0.3 for x = 0.30 in the well- 
ordered magnetic state but also increased with doping, 
i.e., for X = 0.25 Ar/(27r/A) ~ 0.04 at T « 5 K (visible 
in Fig. la). The increase of the damping rate with dop- 
ing indicates that Na doping indirectly causes disorder in 
the FeAs layers, albeit much weaker than Co doping that 
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FIG. 3. (a) We use TP fiSK to estimate the paramagnetic 
volume fraction of Bao.7oNao.3oFe2As2 for Tc2 < T < Tn to 
scale the ZF precession amplitude and obtain the magnetic 
volume fraction shown as circles in (b), the phase diagram 
of the Bao.7oNao.3oFe2As2 sample: Below Tjv = 40.0(5) K, 
77(2)% of the sample volume developed long range antifer- 
romagnetic order. Below Tc2 ~ 29.0(6) K, the same volume 
also became superconducting, which was accompanied by a 
reduction of the ordered magnetic moment by ~ 65%. The 
remaining 23(2)% sample volume stayed paramagnetic and 
became superconducting below Tci = 34.5(13) K. 



causes overdamped muon spin precession close to opti- 
mal doping [13]. Because Na is introduced between the 
FeAs layers, whereas Co resides directly in the magnetic 
FeAs layer, such behavior be expected. We summarize 
the results of TF- and ZF-/.jSR for x = 0.30 in a phase 
diagram shown in Fig. 3b. 

The temperature dependence of the precession fre- 
quencies /yi(T), fsiT) shown in Figs, lb, c are identi- 
cal and the temperature dependencies of the precession 
amplitudes, shown in Fig. 3b for x = 0.30 (scaled to the 
magnetic volume fraction), arc constant below w 0.95T/V 
for both magnetic compounds. This observation indi- 
cates a magnetic structure that remains unaffected by 
superconductivity. In particular, it indicates that the re- 
duction of the precession frequencies for x = 0.30 below 
Tc2 is solely due to the reduction of the ordered moment 
and not due to a reduction of the magnetic volume. With 
a local probe like /iSR this distinction is possible, whereas 
scattering techniques arc only sensitive to the product of 
both quantities. 

The magnetic order for x — 0.25, and 0.30 prohibits 
the study of the penetration depth by /iSR. For x = 0.40 
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FIG. 4. The in-plane penetration depth \^^{T) (proportional 
to the superfluid density) and the (diamagnetic) average in- 
ternal field of Bao.6oNao.4oFe2As2 determined by TF /iSR. 



we measure the probability distribution P{B) of local 
magnetic fields in the vortex state by TF-/.tSR with 
I-IqH = 15 mT parallel to the crystallographic c-axis. 
Fig. 4 shows \~^{T), which is proportional to the stan- 
dard deviation of P{B), and the average local magnetic 
field [31]. We find that X^^{T) is best described by two 
independent s-wave gaps without nodes A(0) = 8.3(5), 
3.0(4) meV, = 33.5(1) K, the weight of the large 
gap is w = 0.66(6), and 1 — cj for the small gap and 
Aa6(0) ~ 262 nm (see supplemental material for de- 
tails). Vortex lattice disorder would artificially reduce 
\ab (broaden P{B)), therefore Xab{0) ~ 272 nm is a lower 
limit for the in-plane penetration depth. These values 
are in good agreement with A = 10.5, 3 meV determined 
ARPES on a single crystal from the same batch with 
the same composition [21]. Specific heat measurements 
[23] on a crystal with 35% Na doping and Tc = 29.5 K 
revealed A = 5.27, 2.68 meV, which indicates that the 
small gap increases linearly with Tc {2A/ [ksTc) = 2.1 for 
both compositions), whereas the large gap grows non- 
linear with Tc {2A/{kBTc) = 4.2,5.8(4) for x = 0.35, 
0.40, respectively). 

In summary, we have conducted muon-spin relaxation 
experiments on Bai_a:Naa;Fe2As2, with x = 0.25, 0.30, 
and 0.40. The Na doping reduces the magnetic ordering 
temperature to Tn = 123(1), and 40.0(5) K for a; = 0.25, 
and 0.30, respectively, but the magnetic transitions re- 
mains sharp and the magnetic order well defined. For 
X = 0.40, magnetic order and slow fluctuations were ab- 
sent down to T « 2 K, as opposed to optimally doped 
BaFe2-a;Co2:As2 [13]. All three specimens are bulk super- 
conductors with critical temperatures Tc = 9, 34.5(13), 
and 33.5(2) K for a; = 0.25, 0.30, and 0.40, respec- 
tively (Tc for X = 0.25 from Ref. [21]). For x = 0.25, 
and 0.30 we presented unambiguous evidence for micro- 
scopic coexistence of the two phases. The reduction of 
the magnetic order parameter by superconductivity was 
strongly dependent on doping: it was too weak to detect 
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for X = 0.25, but amounted to w 65% for x = 0.30. For 
underdopcd Bai_2:Ka;Fe2As2; on the other hand, Wiesen- 
mayer et al. [17] reported a much smaUer reduction of 
the ordered magnetic moment below that did not in- 
crease significantly with doping. Within our itinerant 
model this can be accounted for by a much weaker su- 
perconducting order parameter relative to the magnetic 
one, i.e., better nesting of the electron and hole bands. 
On the other hand, the non-monotonic doping depen- 
dence points to disorder effects which could suppress 
magnetism and superconductivity differently [32]. The 
presence of a second, purely superconducting minority 
phase for x = 0.30 suggests chemical inhomogeneity — 
nevertheless, microscopic coexistence occurred in the ma- 
jority of the sample volume (77(2)%), suggesting that 
this phenomenon is more robust against disorder and 
doping in this system than in Bai_2:Ka;Fe2As2. The sym- 
metry and size of the superconducting gaps, as well as 
Tc for X = 0.40 is similar to those of Bai_2:Ka;Fe2As2 
[33]. Considering the evidence for microscopic coexis- 
tence presented by Wiesenmayer et al. [17] in under- 
doped Bai_a;Ka;Fe2As2 polycrystals and in this Letter for 
under- and nearly optimally doped Bai_a;Na2:Fe2As2, we 
conclude that this phenomenon is intrinsic to hole doped 
BaFe2As2 and can be described by an itinerant two-band 
model as for BaFe2_a;Co^As2 [26]. 
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Analysis of \-^{T) 



In analogy to Ref. [33] we analyzed y^^{T) assuming 
two independent contributions to the superfluid density, 
arising from two s-wave gaps without nodes with different 
gap sizes Ai(0) and A2(0) but identical Tc- We let the 
gap size independent of Tc, which is the so called a-model 



[34]. We fit the following expression to the data [35, 36] 

Aab(T)-2 

(1) 



D{AiT),T) =2 / 

J/5 



XabiO)-'' HI - D{Ai{T),T)) 
+{l-oj){l~D{A2{T),T))], 
d,f{E) E 
A dE ,/W~K{Tf 



dE, 



where f{E) is the Fermi-Dirac distribution, to the weight- 
ing factor of the gaps, and the temperature dependence 
A(r) given by Carrington and Manzano [35]. 



